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Abstract 

We describe a calculation using quenched lattice QCD of the connected part of the neutron 
electric form factor for momentum transfers in the range 0.3 GeV'^ ~ —cp' ~ 1.0 GeV^. We extract 
the implied charge radius using a Galster parameterization and consider various ways of extrapo- 
lating the neutron charge radius value to the chiral limit. We find that the measured charge radii 
may be reconciled to experiment by standard phenomenology and lowest or next to lowest order 
contributions from chiral perturbation theory. 
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I. INTRODUCTION 



There has been a surge of interest and activity in evaluating electromagnetic form factors 
of nucleons, both on the experimental and theoretical (lattice) side S]. Unfortunately, 
knowledge of the neutron electric form factor, Gg(g^), has lagged behind. The extraction of 
this fundamental quantity is difficult because of the dominance of the magnetic contribution 
in unpolarized measurements. In addition, free neutron targets do not exist making its 
measurement from deuteron targets prone to model dependent systematic errors of the 
order of ~ 30% ^. However, new types of experiments with a polarized electron beam 
and deuteron or ^He targets now lead to measurements of G^{q^) with greatly improved 
accuracy 3,0]. 

There are reasons to believe that this form factor is the most interesting and revealing of 
the four nucleon electromagnetic form factors. The contribution of the various flavored sea 
and cloud quark components will no longer be hidden behind a huge valence contribution, 
as is the case for the other electromagnetic form factors, and the sizes of the non-valence 
components could hold some surprises. Lattice QCD was invented to reliably sort out such 
subtle effects, and it will be a major accomplishment to predict /explain these properties. 

The calculation described here is of the "connected" part of G" (g^). The electromagnetic 
current is a sum of the charge- weighted flavor currents, each of which is a color singlet. 
Thus there are also quark self-contraction graphs, usually called "disconnected" graphs, 
involving the various flavor currents which can also contribute to this quantity. (These self- 
contraction loops are of course connected to the valence or cloud quarks by gluons.) We do 
not attempt to calculate these graphs, but will comment on the signiflcance of our results 
for the disconnected contribution in Section V. 

Previous lattice calculations of nucleon electromagnetic form factors using Wilson ^ 0, 
Q| and improved fermions |2] have had good signals for the proton electric and magnetic 
form factors, Ge(q'^) and G'f^(g^), as well as the neutron magnetic form factor, C^i^q^), but 
noisy results for the neutron electric form factor, G^^{q^). Similar to experiment, the lattice 
extraction of this form factor is the most difficult. If lattice studies are to keep pace with the 
improved experiments referenced above, we must move beyond the qualitative stage in the 
calculation of this physical quantity. A new, more precise, lattice calculation of G" (g^) is 
described here. The lattice used is larger (20^ x 32), more configurations (100) are employed 
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than in Ref. [5|, and more sophisticated analysis methods are used. 

Although we will calculate the electric form factor at several nonzero momentum trans- 
fers, we prefer at this point to concentrate on the implied charge radius, a zero momentum 
quantity which should allow a more reliable chiral extrapolation. We will see that the neutron 
charge radius is a useful laboratory in which the physics of the quenched approximation can 
be studied, and its shortcomings quantified. One of the effects of quenching on our results 
comes through the double pole "hairpin" graphs. We will consider these graphs separately 
and provisionally conclude that these do not have a significant effect on our extrapolations. 
We think it extremely likely that the quenched approximation gives a good representation 
of the connected part of G^{q'^). The quenched, connected simulations quoted above have 
yielded phenomenologically acceptable chirally extrapolated electromagnetic form factors 
for both mesons and baryons. The quenched approximation has also been generally suc- 
cessful in evaluating other quantities such as, for example, scalar and axial nucleon form 
factors The validity of this approximation, however, can not be assured and subtle 
dynamical/disconnected effects may remain. We will suggest further investigations which 
should help form a more complete physical picture in Section V. 



II. NUMERICAL DETAILS 



The periodic gauge configurations used in this study were generated from the standard 
one-plaquette action at /? = 6.0 on 20^ x 30 lattices. These lattices were slightly enlarged to 
20^ X 32 by copying each time edge to the opposite side because of an inverter requirement 
that at least one of the dimensions be an integral multiple of 8 0]. The Cabibbo-Marinari 
pseudo heatbath was used and configurations were separated by at least 1000 sweeps. No 
significant auto correlations were observed in the real or imaginary parts of the nucleon 
correlators. 100 configurations were used in the charge radius analysis. 

We used k = 0.150,0.152,0.153, and 0.154 in the standard Wilson quark action. The 
pion masses used are 



m^a = 0.578(2), 0.47721^ |10|, 0.4237(8) fill, 0.364(1) 



0, 0.4237(8) ju']. 
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respectively. The nucleon masses, previously reported in Ref. jl^ on 2000 configurations, 
are 
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mjvci = 0.997(7), 0.869(2), 0.799(2), 0.728(3) 

(The pion and nucleoli masses at k = 0.150 were measured here on t=15 to 18 single 
exponential fits, where t=0 is the time origin of the quark propagators.) The lattice spacing 
is taken to be the same as in Ref. Q, 

a = 0.1011(7) fm 

obtained in Ref. from the physical string tension. The authors of Ref. j3| used the 
physical nucleon mass instead to arrive at a = 0.115(6)fm. 

The electric form factors were measured from point source neutron two and three point 
functions with the ratio given in Ref. ^J, with zero and non-zero momentum two point 
functions evaluated at timestep t' = 9: 

r ( ( V'"^ / g«J4n(t2,ti;0,-g,r4) \ / Gn„(t^;0,r4) \^ , . 

^^^^ \E + m) \ G„„(t2;0,r4) l\Gnn{t'-^q,T,)j ■ 
The time ordered two point function, using the standard neutron interpolation field, 

X^{^) = ^'^'^'^'ipa^^{x)ipjf^''{x){C)/3^tlj^^'^{x), where C = C'j^ and C is the charge conjugation 

matrix, is (understood a, a' sums) 

G'„„(t;p,r) = 5]e-^^^--T,,,(vac|T(xS(x)x^,(0)) |vac), (2) 

X 

I r 



where r4 = ^ 



in the 4x4 Dirac space. The three point function we need, which 







uses the lattice conserved vector current, J^{x), is (g = p — p' 



V 



Gr,j^^{h,h-p,p',T) = -zY,e-'P^-^'e'^^^^^^ (3) 

X2,Xl 

The time position of the conserved charge density operator, J4{xi), which is extended one 
lattice link in the time direction, is associated with the midpoint of that link, which means 
half-integer values for ti. 



III. RESULTS 



In Figs. 1(a) and (b) we show an example of the raw data fits, in this case for k, = 0.152 
and the lowest momentum transfer, where the time axis represents ti in Eq.(0). In order 
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to establish reliable source positions for maintaining a good signal, we originally did the 
calculations on 50 configurations with point nucleon sources at time steps 6 and 27. This 
is the At = 21 data referred to in Fig. 1(b). The (b) part of the figure shows that a good 
plateau is already forming at time step 7, although error bars are large. This allowed us to 
move the final source position to time location 21, giving At = 15 between sources, and three 
point correlation data at time steps 7-10 relative to the origin were fit. The lowest three 
momentum transfers had similar plateaus and form factor values were extracted for the other 
K, values. Fig. 2 illustrates the form factors found at k = 0.154 and 0.153; numerical values 
and error bars are given in Table I. The error bars on the form factors only appear relatively 
large because of the small central values. Comparing results at k = 0.154, the error bars 
here are a factor of approximately 5 times smaller than the similar calculation in Ref. P|. 
(They are also significantly smaller than error bars given in Ref. j^, although individual k 
results are not given.) With the adopted scale for a, our smallest four- momentum transfer 

n 

range on these measurements is about 0.3GeV^, about half that used in Ref. p. 

Before we go on to the charge radius fits, let us explain some of the philosophy of our 
fits of the electric form factor data. A well known phenomenological form for G^{q^) is the 
Galster parameterization 

GT\q') = -YT^G^il')^ (4) 
where r = —q'^/4:m%, where is the nucleon mass and 

G^(g2) = (5) 

is the dipole form factor using the dipole mass, rriu. The most often used value of the 
parameter p is 5.6 and fi ~ —1.91 is the neutron magnetic moment. The neutron squared 
charge radius, given by 

-n = 6^^|,..o, (6) 

imphed by the Galster form is = 3/i/2m^ = —3.24 GeV~^ which is only about 9% 
different from the experimental value, = —2.982(6) GeV~^ [l5|. This form is used to fit 
the form factor data given in Table I. At each value of n we have three values of G'g(g^). 
There are also three parameters in the Galster form: p, /i, and rriD. It is essential to have 
at least one degree of freedom in the fits so that error bars and goodness of fits can be 
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defined. Dipole masses for Wilson fermion fits of tlie proton electric form factor were given 
in Table VII of Ref. 0] at k = 0.154 and 0.152. These values are adopted as input and listed 
in GeV in Table II, along with the nucleon masses given above. In addition, the values 
at K = 0.153 and 0.150 were interpolated from the n = 0.154,0.152,0.150 dipole masses 
given in Ref. P]. These were plotted as a function of dimensionless quark mass, defined 
as ln(4Kc/3K — 3), where we used the central value in Kc = 0.157096(28)19^ Q|. A linear 
plus constant fit produced an excellent interpolation with = 0.025. The final Galster 
parameters for the four fits are listed in Table II and are shown in Fig. 2. Again, the error 
bars on the extracted -r^ appear large, but even at the largest k the error bar is about 6 
times smaller than the experimental value for this quantity. The error bars could have been 
made significantly smaller by fixing the p value in the Galster fits, but we choose not to do 
this. 

In our chiral fits we are extrapolating the neutron charge radius using the formulas 
in the Appendix from heavy baryon chiral perturbation theory and parameters in Table 
III across our four values to the physical charged pion mass at 139 MeV. In doing 
these extrapolations, it is necessary to adopt values of /^r and A (the octet-decuplet mass 
difference). Our point of view is that we are extrapolating down toward the chiral limit from 
our lowest pion mass at k = 0.154, so our values of and A are measured at this value Q|. 
With our value of a-\ the result of Table 4 of Ref. jlj gives A = 158 MeV, and Table VIII of 
Ref. [l^ gives f-,, = 120 MeV at k = 0.154 when the value adopted here is supplied. (It 
is also consistent with the interpolated f^^ value at k = 0.154 from Table XX of Ref. 10l|.) 
We also adopt the tree level QCD values D = 0.8 18], the axial vector octet-octet coupling, 
and \C\ = 1.4 the axial vector octet-decuplet coupling. These last two quantities are 
fixed from phenomenology at the chiral limit, so we can not claim complete consistency j^]. 
Of course, other extrapolation schemes can be used, and the first reference in j|2| for example 
assumes some parameter values at the chiral limit and extrapolates upward in toward 
their measured (isovector) charge radii and magnetic moments using other parameters. 

Fig. 3 shows two extrapolations of the -r^ values as a function of using these param- 
eters. One curve in Fig. 3 assumes only quenched octet intermediate states, labelled "O" in 
Table IV. (Ci, C2 give the analytic contribution to -r^ and are defined in the Appendix.) 
Perhaps surprisingly this produces an excellent fit of the data. Although the charge radii 
values are initially small, they extrapolate to a physical charge radius which is consistent 
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within errors with the measured experimental result because of the chiral logs. The other 
curve in Fig. 3, the "0+D" result from Table IV, gives the result for adding the leading 
octet and decuplet contributions. As one can see this lowers the extrapolated -r^, which is 
still consistent with experiment. The error bars on the charge radius values at the chiral 
limit come almost completely from the uncertainty in the Ci values from Table IV, since 
the physical pion mass is essentially at the chiral limit. The uncertainties in Ci and C2 
(given by ±0.43 GeV"^ and ±0.35 GeV""^, respectively) are independent of the fit (as is the 
uncerta inty in the extrapolated -r^, as listed in the Table) since in the Marquardt algorithm 
of Ref. [2l| the uncertainties are determined by either analytical or numerical derivatives of 
the functional form with respect to the adjustable parameters. Since these are just constant 
plus linear fits in m^, the fixed nonanalytic part does not affect the uncertainties in Ci, C2 
or the extrapolated -r^. 

Fig. 4 contrasts the "0±D" fit found in Fig. 3 with two changed forms which also involve 
both the octet and decuplet. One of these is the "cutoff' form suggested in Ref. j3| in 
the context of parton distribution functions, labelled as "(0±D)(c)" in Table IV. The other 
form is the higher order (in l/rriheavy) nonanalytic calculation of Ref. j^], adapted to this 
quenched situation and labelled "(0±D)(h)" in Table IV. The figure illustrates that these 
forms can lower the extrapolated -r^ and can have a significant impact (~ 10 — 30%) upon 
the chiral limit, which we think is a reasonable estimate of the systematic extrapolation 
error involved. The "(0±D)(h)" extrapolation raises the value from the "(0±D)" fit and 
essentially hits the experimental value exactly with a lowered Xr- would be wrong to 
make too much of this agreement, but clearly the fits we are making are consistent with the 
lattice data and the experimental value for the neutron squared charge radius. 



IV. HAIRPIN CONSIDERATIONS 



We have considered the quenched QCD contributions to the neutron and proton charge 
radius from the so-called hairpin graphs using the methods of 2M- The leading nonanalytic 
parts of these diagrams are proportional to 51n(m^/yU^) where 5 = ml/{487i^ f^) (for m„ = 
ma) and rriQ is the usual hairpin "mass" parameter. These diag rams are the direct analog 
of the magnetic hairpin graphs considered recently in Ref. 25[. (See Figs. 2(a) and (b) 



of this reference for the relevant diagrams. The a^k interaction should not contribute to 
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a leading log for magnetic moments or charge radii because of the two extra derivatives 
in the coupling.) Remarkably, we find that the contributions from these graphs vanish 
identically for the proton and neutron. There are three sets of such graphs plus wavefunction 
renormalization. The first set consists of a hairpin correction of the tree-level octet charge 
radius vertex. When this is combined with the wavefunction renormalization, these two 
contributions cancel identically, not only for the proton and neutron, but also for the other 



octet baryons 



26|. In addition, there is a contribution at this order to the charge radius from 



the electric quadrupole operator [27[. It connects octet and decuplet baryon lines similarly 
to the transition magnetic moment operator in Ref. 2^. It vanishes for the proton and 
neutron (although not for the other octet baryons) for the same reason that the magnetic 
transition operator does; namely, that the flavor-charge coefficient vanishes. In addition, 
there are one-loop hairpin contributions from the tree-level decuplet vertex, which also 
includes a wavefunction renormalization part. These do not cancel, but vanish separately 
for the proton and neutron, again because of the flavor-charge factor 3]. 

Unfortunately, there are also two-loop graphs which contribute at the same order as the 
one loop graphs within chiral perturbation theory. This can be understood simply from 
the physical dimensions of the charge radius, which goes like (mass)"^. To lowest order in 
the hairpin interaction the modification of the tree-level charge radius is of order 5(r^)tree- 
However, a term of the same order with the usual chiral perturbation theory factor l/(47r/^)^ 
replacing the factor of {r'^)tree niay also be generated from a two-loop graph involving a 
primitive electromagnetic interaction. Two-loop graphs have been considered recently in 
the context of mass corrections in Ref. [2^. These graphs are clearly related to the graphs 
we need since one of the pion loops may be changed to a hairpin loop just by differentiation 
with respect to the pion mass and changing the flavor factor. Although we have stopped 
short of their evaluation we can conclude from power counting that these charge radius 
graphs are logarithmically divergent in the loop momentum and can in principle contribute 
to the chiral log charge radius coefficient. 

An estimate shows that the contribution of any such hairpin chiral logs in the proton and 
neutron charge radii should be small. Recently, there have been estimates of the value of 6, 
the hairpin coefficient, from lattice simulations 3^. For concreteness we will take S = 0.2. 
As an example of the expected magnitude of the effect of the delta interaction on the chiral 
log coefficient, let us take the magnitude of the ratio of the wavefunction renormalization 
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contribution (canceled by the modified tree-level piece) to the quenched chiral coefficient for 
the proton or neutron: 

R ^ {6{rl),UD - n')/{j^,fD')- (7) 

Our fits below (see the Ci coefficients in Table IV) give (r^)free < 1 GeV^, so that for D = 0.8, 
F = 0.5 [18], we obtain at most R ~ 0.05. If we replace the [rDtree by another factor of 
1 / (AirfTj)^, which would happen at two loops, we get the same estimate. Of course, the flavor- 
charge coefficient here is specific to this interaction. However, barring an exceptionally large 
coefficient it is likely that the effect of hairpins on the chiral extrapolations will be negligible. 



V. SIGNIFICANCE OF RESULTS 



We have seen that one can come remarkably close to understanding the value of the 
neutron charge radius using quenched lattice QCD. Our calculations use the Wilson action 
and therefore are limited to rather large pion masses. Although standard values for the D 
and C coefficients were used in the analysis, there is nothing to force the quenched values 
of these coefficients to be the same as in the full theory. The values we used should be 
considered merely typical, based on standard phenomenology. However, our results for -r^ 
are consistent with the expected rise from the quenched chiral logs using these coefficients, 
and the ffis to the data are excellent. Of course lattice calculations, both quenched and 
dynamical, probing further into the chiral limit would be helpful to verify this scenario. 

We find that the higher order nonanalytic contributions in Ref. j^] or the possible "cut- 
off' form suggested in Ref. j2^] can make a small but significant impact on the chiral extrap- 
olations. According to the argument in Section IV, we do not expect hairpin contributions 
to make a significant contribution to the chiral log terms. However, we can not rule out 
the possibility of large coefficients entering and this remains a slim possibility. An exact 
evaluation of this quantity would be very welcome. 

Of course, the disconnected part, which consists of up, down and strange quark loops, 
should be added to our results before we compare to experiment. Ref. 13] gives results for 
combined u, d, and s quarks and finds an increase of perhaps 0.5-1.0 GeV~^ ~ .019 — 039 fm^ 
in — r^. (We are using the smallest — result from Fig.5(b) of this paper and the Galster 
parameterization.) The question naturally arises as to whether our results are compatible 
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with such values. The fact that one finds good agreement with experiment using standard 
phenomenology and only the connected part of the amplitude argues that the disconnected 
contribution must be relatively small, again assuming the hairpin contribution is minor. An 
increase by the amount suggested in Ref. |^ could perhaps be accommodated by our data, 
but likely at the expense of the excellent fits found here. The present results illustrate the 
importance of establishing a reliable set of extrapolation parameters, which can only be done 
by extending lattice calculations for a variety of quantities closer to the chiral limit. For 
now, we conclude from our lattice study that the connected part of the quenched amplitude 
is capable of explaining the bulk of the neutron electric charge radius. 

Note added: An independent discussion of charge radii in quenched and partially 
quenched ChPT has just appeared 0]. 
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Appendix A: Quenched Chiral Charge Radii Expressions 



We use the ~ 93MeV normahzation. For completeness we give the lowest order chiral 
expressions for both octet baryon charge radii as well as magnetic moments. The charge- 

1 '^N 1 



flavor factors, fei'f'', are listed in Table III. 



-Aln(^)+2v'm|-A^cos-(;^) , > A 
-A In ( ^ j - 2^A2 - m\ ln( — — ^) , A > 



(*^'^).-(^ E/?'(-|^4™x,A,.)-f) (12) 



1„ f ^ ) + cos-' f ^1 , > A 

,2 



F,(mx,A,/i) = <; )^;( , >x^^ (13) 
In ^ + , 1^ , ln( — ^ , A > mx 

V '' pi. a/a ™x ^ mx ^ ' 

Note the sign typo in Eq.(6) of Ref. 23[ in the A > mx branch of Fr{mx, A, /i). Although 
these quenched expressions do not appear in the literature, it is clear that they are implied 
by the results in Ref. js^ for the full QCD charge radii and the results of Ref. for the 
quenched magnetic moments since one can simply take the ratio of the quenched to full QCD 



coefficients in Table III of Ref. |25|, and apply them to the full QCD charge radii results in 
js^. (One must also drop the D, F independent terms from Table 1 of |32|.) We have also 
independently verified the above results for charge radii and magnetic moments. We have 
left analytic parts in the expressions so that others may independently verify these results 
using dimensional regularization. 

The functional form of the nonanalytic expressions used in the chiral extrapolations of 
the neutron — are given below. 
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Octet only: 

Octet plus decuplet: 

Octet plus decuplet cutoff form: 

1 fin /A2_|_rn2\ cr ^ 

(O + D)(c) = O + D + ^^-^ I -D' In [-^r^j " f'Fri^ml + A^ A, /.)| . (16) 
Octet plus decuplet with higher order (~ 1/heavy^) correction: 
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(Prim^, A, /i) - j^Gim^, A, /x)^ | 



where (Ref. |23l |) 



A In (< ) + ^^=% cos-i (^) , > A 

To all such forms for — are added the analytic terms Ci + C2m^ where Ci, C2 are constants. 
Note that there can be additional terms depending on A in the "cutoff' form, Eq.(|16p. but 
we include only the modifications to the nonanalytic terms in Eq. ()15|l in the same spirit as 
Ref. y. 
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TABLE I: Numerical results on neutron form factors at four values of k and various values of the 
four momentum transfer —q^. 



K 


-q2 (GeV2) 


GJ-q^) 




u.ooo 




0.154 


0.689 


0.0134(43) 




0.996 


0.0640(52) 




u.ou± 


U.UUOOl Z(± I 


0.153 


0.699 


0.0099(29) 




1.02 


0.0054(36) 




0.364 


0.0067(15) 


0.152 


0.707 


0.0076(21) 




1.03 


0.0046(27) 




0.366 


0.0042(8) 


0.150 


0.717 


0.0049(13) 




1.05 


0.0034(17) 



TABLE II: Galster fit parameters. The values of p and — are results of the fits, the values of 
rriN and itid are input, is the chi-squared of the fit. 



K. 


mo (GeV) 


mjv (GeV) 


P 


-rl (GeV-2) 


x' 


0.154 


1.16 


1.42 


23.(36) 


0.65(51) 


0.62 


0.153 


1.24 


1.5G 


30.(12) 


0.18(31) 


0.50 


0.152 


1.32 


1.69 


33.(42) 


0.34(20) 


0.41 


0.150 


1.47 


1.94 


38.(39) 


0.18(7) 


0.27 
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TABLE III: Charge-flavor factors bj^ of baryon and Goldstone boson X of the quenched non- 
hairpin 1-loop diagrams. The symbols B and T denote the octet and decuplet baryons respectively. 



N 








6f 


P 


16 ^2 
3 





1 /72 

2 ^ 





n 


_16 ^2 
3 ^ 





1 /72 

2 ^ 





E+ 





16 £,2 
3 ^ 





1 n2 

2 ^ 


SO 





3 ^ 





1 p2 
4 ^ 
















A 





3 ^ 





1 n2 
4 ^ 







_16 £)2 
3 ^ 





1 /72 

2 ^ 

















TABLE IV: Final fit parameters on some chiral extrapolations of — r^. The parameters in the fits 
are: A = 158 MeV, = 120 MeV, D = 0.8, \C\ = 1.4, /x = 1 GeV, A = 550 MeV. "O" denotes 
octet only; "0+D" denotes octet plus decuplet; "(0+D)(c)" means the cutoff form, "(0+D)(h)" 
means the higher order octet plus decuplet form. Xr = X^/^df is the reduced chi-squared of the 
fit, where the number of degrees of freedom, Ndf = 2. 



fit 


Ci(GeV-2) 


C2(GeV-4) 


-r'n (GeV-2) 







-0.230 


0.501 


3.49(43) 


0.0008 


0+D 


0.060 


-0.093 


2.49(43) 


0.010 


(0+D)(c) 


0.580 


-0.373 


2.22(43) 


0.022 


(0+D)(h) 


0.321 


0.252 


2.90(43) 


0.003 



16 



0.04 




tl 



(b) 



O (5 B 



o n n 



B e o 



1 2 
tl 



FIG. 1: (a) k = 0.152 electric form factor data for p = vr/lO on the source separation At = 15 
lattice for 100 configurations, (b) Similar graph on the source separation At = 21 lattice for 50 
configurations. The fit value on the At = 15 data is shown superimposed on the At = 21 data in 
the (b) part. 
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FIG. 3: Two chiral extrapolations of the neutron charge radius data (in GeV~^) as a function 
of pion mass squared (in GeV^). The solid line is the pure quenched octet ("O" in Table IV) 
extrapolation, the dashed line gives the octet plus deeuplet ("0+D"). The vertical dashed line 
represents the physical pion mass squared and the burst on this line represents the experimental 
charge radius. The lattice data are represented by open circles and are given as a function of 
squared pion mass in GeV^. Error bars are given on the extrapolated values (triangle - octet, solid 
circle - octet plus deeuplet) at the physical pion mass squared. 
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FIG. 4: Three chiral extrapolations of the neutron charge radius data (in GeV~^) as a function 
of pion mass squared (in GeV^). Prom top to bottom (shortest dashed Une to longest). All fits 
take into account both the octet and dccuplct. The top line gives the result of the higher order 
fit ("(0+D)(h)" in Table IV), the middle is the (lowest order) octet plus decuplet ("0+D") also 
shown in Fig.3, and the lowest dashed line is the cutoff form ("(0+D)(c)"). The solid square, 
circle and triangle give the chiral values of these three fits, respectively. The solid circle has been 
moved to the left for clarity of observations of the various error bars. See Fig. 3 for meanings of 
the other symbols. 
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